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We present fundamental constraints required for a consistent linear response theory of fermionic 
(~| ■ superfluids and address temperatures both above and below the transition temperature Tc- We 

' emphasize two independent constraints, one associated with gauge invariance (and the related Ward 

identity) and another associated with the compressibility sum rule, both of which are satisfied in 
strict BCS theory. Ifowever, we point out that it is the rare many body theory which satisfies 
both of these. Indeed, well studied quantum Hall systems and random-phase approximations to the 
electron gas are found to have difficulties with meeting these constraints. We summarize two distinct 
theoretical approaches which are, however, demonstrably compatible with gauge invariance and the 
compressibility sum rule. The first of these involves an extension of BCS theory to a mean field 
^j^' description of the BCS-Bose Einstein condensation crossover. The second is the simplest Nozieres 

Schmitt- Rink (NSR) treatment of pairing correlations in the normal state. As a point of comparison 
we focus on the compressibility n of each and contrast the predictions above Tc- We note here that 
^ ' despite the compliance with sum rules, this NSR based scheme leads to an unphysical divergence 

, in K at the transition. Because of the delicacy of the various consistency requirements, the results 

^H. of this paper suggest that avoiding this divergence may repair one problem while at the same time 

introducing others. 

'a " 

PACS numbers: 03.75.Ss,74.20.Fg,67.85.-d 

C . I. INTRODUCTION: THE CHALLENGES OF ARRIVING AT CONSISTENT LINEAR RESPONSE 

O ■ THEORIES 

With recent progress in studies of ultracold Fermi superfluids undergoing BCS-Bose Einstein condensation (BEC) 
crossover has come a focus on generalized spin and charge susceptibilities (see, e.g., Ref. [l| for an introduction). 
^ These in turn relate to the compressibility and spin susceptibility in the linear response regime. While strict BCS 
theory is well known to lead to fully consistent results for linear response (see Ref. 0] for a review), arriving at a 
generalization to address the entire crossover is a major challenge. The difficulty is to construct such theories as to 
be fully compatible with /-sum and compressibility-sum rules which reflect conservation principles. Although there 
have been some successes there are nevertheless important failures, as we will address here. 

In this paper, we discuss these challenges in the context of fermionic superfluids. We present examples of theoretical 
I approaches which are fully consistent with /-sum and compressibility-sum rules. To our knowledge there are two such 
theories beyond strict BCS theory and we discuss both here. The first of these involves a mean-field approach to 
BCS-BEC crossover. Here we avoid the complexities associated with pair fluctuation or pseudogap effects The 
second of these involves an investigation of the normal phase which includes pairing fluctuations at the simplest level. 
In this context we address Nozieres and Schmitt-Rink (NSR) [4] theory. We demonstrate that gauge invariance and 
^ ' sum rules can be made compatible with a linear response scheme. Importantly, however, this linear response theory 
is unphysical in that it predicts a divergent compressibility when approaching the superfluid transition temperature 
from above. This observation should underline the point noted above that a treatment of linear response functions in 
fermionic superfluids (and many-body systems overall) is full of pitfalls. 

Indeed, for a general many-body theory, as summarized in Ref. [Sj], it is difficult to obtain the same result for the 
compressibility via the derivatives of thermodynamic quantities as compared with that via the two-particle correlation 
functions. This issue has been extensively discussed for random-phase approximation (RPA) and RPA-generalizations 
of the electron gas Q as well as in quantum Hall systems ,7]. From this literature, it appears that when the two 
approaches for the compressibility disagree, the more trustworthy scheme Q is the thermodynamical approach. In 
many ways more complex are the response functions of BCS theory which from the very beginning \§, ^ revealed 
difficulties associated with incorporating charge conservation and gauge invariance in the broken-symmetry phase. 
In a related fashion we note that there are considerable discussions about where and when collective mode con- 
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tributions (which in neutral systems are phonon-like) associated with the fluctuations of the order parameter enter 
into the electromagnetic response and related transport of fermionic superfluids. These phononic modes are central 
to the bosonic superfluids. Because correlation functions (and related transport coefficients) are constrained by sum 
rules, one is not free to incorporate these collective modes in an arbitrary fashion. In this paper we discuss in some 
detail how these collective modes appear in a theory respecting conservation laws. We will show that in the process of 
maintaining gauge invariance one must self-consistently calculate the collective mode spectrum. The usual phonons, 
which appear as density waves above the transition, are below strongly entangled with the phase and (in general) 
amplitude modes of the order parameter. 

Indeed, it appears that in the BCS superconductivity literature (including helium-3 as the neutral counterpart), 
phonons do not contribute to transport properties associated with the transverse correlation functions. This includes 
the conductivity, and shear viscosity. A consequence of this body of work based on kinetic theory as well as Kubo-based 
systematic studies approaches [lol . IT]| is that the normal fluid or condensate excitations involve only the fermionic 
quasi-particles. In superfluid helium-3, as well, the normal fluid contributions to thermodynamics reflect only the 
fermionic quasi-particles. 

This brings us to the possibility that the situation is different for strong coupling superfluids in the sense of BCS- 
BEC crossover. Work by our group showed that within a generalized BCS-like description the sum rules are satisfled 
without including collective mode excitations contributing to the normal fluid density [2] . But there may be alternative 
theories where at strong coupling the sound waves are important in this transport. Indeed it has been argued that 
these sound mode contributions are important in the thermodynamics at unitarity, but it should be noted that they 
were similarly invoked in the BCS regime, in a manner which does not appear consistent with theory and experiment 
on superfluid helium-3 [l3| . 

Closely related to the charge and spin susceptibility are the dynamical charge (C) and spin [S) structure factors: 
Ssj Sc- They are formally connected to the charge and spin response functions by the fluctuation-dissipation theorem 
hj. In a proper theory, the conservation laws both for charge and spin yield well known sum rule constraints on 
Ss,Sc- In recent literature there has been an emphasis on these structure factors [Tsl - fTsl . in part because they can be 
directly measured via two photon Bragg experiments and in part, because they are thought to reflect on an important 
parameter for the unitary gases: the so-called "Contact parameter" . 

In this context it has proved convenient to dcflnc 

5±(c.,q) = i[5c(c.,q)±5s(^,q)]. (1) 



It should be noted that in the literature [15| there is a tendency to decompose these spin and charge structure factors 



into separate spin components so that the density (or spin) response is related to correlation functions of the form 

< [Pt ± Pl][Pt ± Pi] > (2) 
Presuming S^^ — and S-f^ = S^-f then one infers 

S{k,uj)^2[Sn{k,Lo) + S^^{k,Lj)] 

with Saa' — J2n < 0\pa{k)\n X n|p^,(fc)|0 > S{u! — En). In this way the difference structure factor S'_(aj,q) is 
frequently associated with density correlations of the form < p-[ Pi >■ 

This association derives from assuming that all diagrams for the spin and charge response are equivalent and given 
by the same combinations of charge density commutators (with only simple sign changes) involving p^{r) ±p^(r). We 
emphasize in this paper that this is specifically not the case below Tc as a result of collective mode effects which only 
couple to the density response function but decouple from the spin response function. Ref. 0] clearly demonstrates 
this difference for strict BCS theory. Moreover, it needs not generally hold when there is a different class of diagrams 
required above Tc in the spin and charge channels to insure the /-sum rules. While these crucial collective mode 
effects are sometimes inadvertently omitted [17l] in analyzing density-density correlation functions, they are essential 
for satisfying the longitudinal /-sum rule. 

Of interest is a claim in Ref. flGj that the static structure factor (which involves an integral over all frequencies) 
at large wavevector measures the so-called Contact parameter. A rather different observation was made by Son and 
Thompson [l8| who showed that it is the high frequency, large wavevector structure factor which is associated with 
the Contact. More precisely. Son and Thompson il8| investigated the relation between the structure factor and the 
Contact interaction noting how delicate this issue is and that "care should be taken not to violate conservation laws" . 
This is the philosophy at the core of the present paper. 
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II. SUPERFLUID LINEAR RESPONSE FORMALISM 



We begin our discussion of linear response theory with the fundamental Hamiltonian for a two-component Fermi 
gas interacting via contact interactions 

H = J d3xV;t(x)(|^ -//)^,(x) -g J d3x^|(x)^|(x)^^(x)7A^(x). (3) 

We assume the interaction is attractive and g is the bare coupling constant. Here we adopt the convention e = c = 
h = 1 and the metric tensor g'^'^ is a diagonal matrix with the elements (—1,1,1,1). 

The goal of linear response theory is to find the full electromagnetic (EM) vertex associated with the EM response 
kernel K^'^{Q). In the presence of a weak externally applied EM field with four- vector potential A'' = (0, A), the 
perturbed four-current density (5 J'^ is given by 

5J^^{Q)=K^^-{Q)A,{Q). (4) 

By introducing the bare EM vertex ^^^{P + Q,P) = {I, 2±i) and full EM vertex F^(P -I- Q, P), the gauge invariant 
EM response kernel can be expressed as 

K^'^iQ) = 2 y F^(P + Q, P)G(P + Qh^iP, P + Q)G{P) + -/i'^^ (5) 

m 

where h'^'^ = —g'^'^{l — g'^^)- Throughout we define Q = ^ (ir2;,q) which is the 4-momentum of the external field 
with being the boson Matsubara frequency, and P = = (icj„,p) is the 4-momentum of the fermion with aj„ 
being the fermion Matsubara frequency. G(P) is the single-particle Green's function. The "bare" Green's function is 
given by Go(P) = (ia;„ - ^p)-i with Cp = ^ - 



A. Central constraints 



Gauge invariance and conservation laws impose an important set of constraints on any linear response theory. The 
full EM vertex F'^ must obey the Ward Identity Q 

qf,T^{P + Q,P) = G'\P + Q)-G-^P) (6) 

and the identity associated with the compressibility sum rule. The former leads to the gauge invariant condition of 
the response kernel q^if "(Q) = which further leads to the conservation of the perturbed current q^5J^ = 0. The 
compressibility sum rule imposes an identity which we call the "Q-limit Ward Identity " [l^ 

.»r«(P + 0,P)l.- = ^^i-^, (7) 

where S(P) is the self energy, and the relation G^^{P) = Gq^{P) — S(P) has been applied. This will guarantee the 
compressibility sum rule dn/d^ = — -ftr'"^(0,q 0) with the compressibility given by k = n~'^{dn/dfi), as will be 
shown shortly. 



B. Linear Response of BCS superfluids 



In BCS theory of fermionic superfiuids, the order parameter is given by 

A(x)-.g(^^(x)^^(x)). (8) 
In the mean-field approximation, the Hamiltonian in the absence of external fields may be written as 

yd3xV't(x)(^|!._^)^^(x)_yd3^(^A(x)^|(x)V](x)+h.c.). (9) 

There are many reviews [l^, [2l| on how to derive K^'^{Q) at the BCS level. Here we set up an approach which we 
refer to as the consistent fluctuation of the order parameter (or CFOP) theory. Importantly, here, in contrast to the 
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approach of Nambu the changes in the phase and amphtude of the order parameter associated with the external 
fields enter as additional components of the perturbation theory. It is useful, however, to cast this CFOP theory in 
the Nambu formulation. We define a± = ^(cti ± 1(72) where tXi is the Pauli matrix, and introduce the Nambu-Gorkov 
spinors 



V'tF 



^l = [^l,i^^.p]. (10) 



In the mean-field BCS approximation, the Hamiltonian (jH) in the presence of electromagnetic (EM) fields can be 
rewritten in the Nambu space as 

H = T. *i.eP^3*p + *l+q( - ^^^^ + - Aq^+ - Alqa_)vl/p, (11) 

p pq 

Here the order parameter is generalized to include fluctuations from its equilibrium value with Aq = 
f/ ^p(^|,(7_vl/p_|_q)^ which will be imposed as a self-consistency condition. When the external EM field is applied, 
the order parameter is perturbed and deviates from its equilibrium value. The order parameter in equilibrium is A, 
which is at q = and can be chosen to be real. Denote the small perturbation of the order parameter as Aq so 
Aq = A + A^. By introducing Aiq = -(A^ + A'!q)/2 and Aaq = -i(A^ - A'_!q)/2, the Hamihonian ^ splits into 
two parts as H = Hq + iJ' with one containing the equilibrium quantities and the other containing the deviation from 
equilibrium. 

- ^ vK^^vI^p, i/'=^vI,^+q(Aiqai+A2qa2+A^qr(p + q,p))vI/p, (12) 
p pq 

p+ — 

where Ep = ^pUa — A(7i is an energy operator and 7^(p + q, p) = (o"3, ^;f-) is the bare EM vertex in the Nambu space. 
Here the perturbation of the order parameter and the EM perturbation are treated on equal footing and this will 



naturally lead to gauge invariance of the linear response theory. The quasi-particle energy is given by Ep — y£,p + A^ . 
The propagator in the Nambu space is 



where 
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G{P) = -. — V+ ■ If ' FiP) = -^pvpi-. ir--—ir) (14) 

lUJn ~ Ep lOJn + Ep \lUJn - Ep lUJn + Ep ^ 

are the single-particle Green's function and anomalous Green's function respectively and Up,Vp = ^(1 ± ■^)- 
By introducing the generalized driving potential and generalized interacting vertex 

*q = (Alq, A2q,A^q)^, S (p -|- q, p) = ((Tl , (72 , (p + Q, p) ) ^ , (15) 

the generalized perturbed current fj is given by 

r/(r, q) - ^(*t,(r)S(p + q, p)vl'p+q(r)) + -S^^H^-'AAt, q). (16) 
p 

where the component 773 = ( J^) denotes the EM current and 771^2 the perturbations of the gap function. This leads 
to a linear response equation in a matrix form 

r/(w,q) = Q{uj,q) ■ *(a;,q) 

/ Qii(cj,q) Qi2(w,q) Q!|'3(a;,q) \ / Ai(w,q) 
= Q2l(w,q) Q22(^,q) %(^,q) A2(c^,q) | . (17) 

VQ^i(..,q) Q^2(^,q) Q^3(^,q) + ^/^^V V^.(^,q) 

The response functions Qij are 

Q,,(r-r^q) = -^(^,[vl/t,(r)S,(p-Kq,p)*p+q(r)vI/t,,^q(r')S,(p^p' + q)V(^')])• (18) 
pp' 
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Using the Wick decomposition {2l|, we obtain 

Q,j{ini,q) = TrT ^ ^ + Q, P)G(P + Q)S, (P, P + g)G(P)) , 



The gap equation leads to 771,2 = — |Ai_2- and using Eq. pT)) . we find 



(19) 



(311Q22 — Q12Q21 



Q11Q22 — Q12Q2I 



(20) 



where Qn = | + Qu and Q22 ^ ^ + Q22- 

The quantity of interest is the EM response kernel K'^'^ , which has the following form in the Nambu space. 



i^^-(Q) =Tr^(f^(P + g,P)G(P + g)r(^,^ + Q)G(P)) 



After substituting Eq. (|20p into our linear response expression we find 



(21) 



(22) 



from which we obtain 



K'^" = Q^" + SK^"" SK^" = Q11Q32Q23 + Q22Q31Q13 - Q12Q31Q23 ~ Q21Q32Q13 



Q11Q22 — Q12Q21 
3 order parameter ; 

From the expression of K^'^ , the full EM vertex r''(P + Q, P) in the Nambu space can be determined. We define 



(23) 



Here Q33 ~ Q33 + j^^^'^ ■ Hence the effects of fluctuations of the order parameter are included in the response kernel. 



H^ 



Q31 


Q21 






Q32 


Q22 


Qii 


Q12 


Q21 


Q22 



Q32 


Q12 


Q31 


Qii 


Qii 


Q12 


Q21 


Q22 



(24) 



From Eq. (|23p . Rf^" can be expressed as 

if^^(g) Tr^([7'^(P + Q,P)-aiH5'(g)-a2H^(Q)]G(P + Q)r(^,^ + Q)G(P)) 
p 

where we have used Eq. (|T9)) to arrive at Q^!^, Qi^ and g23. One can further identify 

f ^(P + g, P) = riP + Q,P)- <yiK{Q) - '^2H^(g). 



(25) 



(26) 



Our ultimate goal will be to demonstrate consistency with the various Ward identities within this BCS formulation. 
In order to proceed we rewrite Eq. (j23p in the form of Eq. ^ . For this purpose it is convenient to define 

H^(g) = -H5'(g) + zH^(g), n^(g) = -H5'(g) -*H^(g), (27) 

With straightforward algebraic manipulations one arrives at 

K^'^iQ) = 2J2 b'^iP + 0, p)GiP + QWiP, p + g)G(p) + n^(g)p(p + Q)Y{p, p + g)G(p) 
p 

+ n^(g)G(P + Q)Y{P, P + Q)F{P) - Y{-P, -P - Q)F{P + Q)Y{P, P + Q)F{P)\ + -h^'r (28) 



Similar expressions for the density response function has also been obtained using a kinetic-theory approach |22| . 
Substituting the expression F{P) — AGq{—P)G{P) into Eg. (1^5)) . and then comparing with Eq.®, we arrive at an 
expression for the full EM vertex 



r'^(p + g,p) = 7^(p + g,p) + An^(g)Go(-p-g) + An^(g)Go(-P) 
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-A2Go(-P)7^(-P, - Q)Go(-P - Q), (29) 

While rather complex, this represents an important result. The second and third terms correspond to the contributions 
associated with collective- mode effects [H , while the fourth term can be identified with the so-called Maki-Thompson 
diagram 0]. 



We remark that in the standard calculation of the Meissner effect and superfluid density [21[ the collective-mode 
contribution is not included. This is because the current-current correlation functions are evaluated and they corre- 
spond to the transverse components of Eq. (|28l) . One can show that the collective-mode contribution cancel in the 
limit = 0, q ^ '^l so the collective-mode effect may be ignored. In contrast, we will show in the following that 
the collective-mode effect for longitudinal response functions is crucial in restoring gauge invariance. 



III. VERIFICATION OF SELF CONSISTENCY IN LINEAR RESPONSE 

In this section we demonstrate that the Ward identities and the Q-limit Ward identity (associated with the com- 
pressibility sum rule) are consistently satisfied at this mean field level. We first present arguments to show that the 
vertex function in Eq. obeys the Ward Identity Contracting both sides of Eq.f^H) with q^, we have 

g^r^(P + Q, P) = G^\P + Q)- G^\P) - 2S(P + Q) + 2E(P) - ^^I-t^^{G^\-P) _ G^\-P - Q)) 

= G'\P + Q)-G'\P), (30) 

which is the desired result. Here we have used the fact that S(P) = — A^G'o(— P) for BCS superfiuids [23|. It can 
be proved analytically that gauge invariance implies that the density-density response function always satisfies the 
/-sum rule (for details, see Ref.[2|) 



I dujuJXpp{uJ,q) ^ n—. 
Here we define the density-density correlation function 



(31) 



TT 



We turn now to the Q-limit Ward identity from which the compressibility sum rule can be derived [19j . We note 
that G-i(P) G(7^(P) - S(P), so that 



dn 



= -2^G2(P)(l-^|^p) =-2^r°(P,P)G(P)7°(P,P)G(P) 



9/i 

K''"{u = 0,(1^0), (32) 



p 



where in the last line, the expression in Eq. ([5]) has been applied. This analysis demonstrates that the compressibility 
obtained via thermodynamic arguments relates to properties of two particle correlation functions. 

The more explicit proof of the Q-limit Ward identity ([7]) is briefly outlined here. Since limq_i.o Go(P + Q)\u]=o = 
Gn{P), we evaluate r°(P -I- Q, P) in the hmit w = and q 0: 

limr"(P + g,P)|^=o = 1-2A limH?(Q)|^=oGo(-P)- A2g2(-P). (33) 

q— J^O q— J^O 

Using E(P) = -A2Go(-P), the right hand side of Eq. is 

1 - - 1 + 2a|^G„(-P) - A2G^(-P), (34) 

where the identity 9^Go(-P) = -Gl{-P)d^,G^^ {-P) = -Gl{-P) has been applied. Comparing Eqs.(I2Sl) and 
one can see that the Q-limit Ward identity holds for BCS theory only when 

(9A 

— = -limH;(Q)U=o. (35) 
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The left hand side can be evaluated by differentiating both sides of the gap equation 1 = 3- J2p ^(P) with respect to 



~da ~ sr A ( l~2 f{E^) r, df(E^] 

Ei\ Bp ^ ^ 

By using the expressions of the response functions given in the Appendix, one can show that 



(36) 



gii(0,q^O) q^o 



lim n?(Q) . (37) 



j=0 



Thus the Q-limit Ward identity is respected, which then guarantees the compressibility sum rule. 

It is of interest to address why an RPA-based approach usually fails to satisfy the compressibility sum rule [1, Q . In 
the RPA approach, one starts with the expression of the density susceptibility called {dn/dfi)o from thermodynamics 
or from a simpler model. By introducing a summation over a series of bubble diagrams, the RPA approach leads to the 
expression {dn/dfi)B.PA = {dn/dfi)o/[l ~ X{dn/dfi)o], where A is a combination of the coupling constant representing 
bubble diagrams and the combinatoric factors for counting the series f^. For finite values of (9n/9/i)o and A, the 
RPA result always differs from (dn/ 811)0 and this makes it particularly difficult to satisfy the compressibility sum 
rule, which requires a consistent expression for dn/dfi. 



IV. APPLICATION TO BCS-BEC CROSSOVER AT THE MEAN-FIELD LEVEL 



We make the important observation that the arguments presented above for the self consistency of linear response in 
strict BCS theory can be readily extended to treat BCS-BEC crossover theory at the mean field level. This level is, of 
course, not fully adequate because it does not include fluctuations associated with non-condensed pairs. Nevertheless, 
it does provide an example of a fully sum-rule consistent approach to linear response in fermionic superfluids. 

At this mean field level, the number density and the order parameter in equilibrium can be written as a function 
of temperature in terms of the equations 



E[l-|^(l-2/M' --T.'-^- (38) 



Because we contemplate arbitrary chemical potentials away from the BCS regime where /i = Ep, this theory applies 
to the whole BCS-BEC crossover regime. Important here is the introduction of the two body scattering length a via 
a renormalization of the coupling constant (24 . [25j 

1 m 



g Ana ^ 2ep 



where a is the s-wave scattering length and ep = /2m. In the BCS limit, a ^> 0~, A — > and Ep where Ep 
is the Fermi energy. While in the deep BEC regime, a — >■ O"*", A — 00 and -A — s> 0. 

IMI 

One of the best measures of a linear response theory is the calculation of the compressibility 

K = n~'^{dn/d^) 

based on the density correlation functions. This is particularly problematic because of the difficulty of finding the same 
answer as found from thermodynamics. It is of considerable interest, then, to establish the form of the compressibility 
in a theory with full (compressibility) sum rule compatibility. At the mean field level the behavior above the transition 
temperature is that of a free Fermi gas. Below Tc, the compressibility either via thermodynamics or via the two body 
density density response leads to 

|^^-if™(0,q^O) (40) 

where if^o ^ qH + SR^^, 
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Figure 1: Compressibility at the strict mean field level and as function of temperature from BCS to BEC. Tp is the Fermi 
temperature, kq — f;;^ is the compressibility for a noninteracting Fermi gas at T = 0. 



Fig[T] plots the compressibility k as a function of temperature for Xjkpa — —1 (on the BCS side), (the unitary 
point), and 1 (on the BEC side). The figure exhibits an expected thermodynamic signature of a phase transition, 
appearing as a discontinuity in the compressibility at Tc. The discontinuity in n, at can be traced back to the 
appearance of collective-mode term 8K^^ which sets in below Tc and is absent in the normal state. It should be noted 
that at Tc, 8K^^ is finite only when w = 0. When uj approaches but does not equal 0, we have bK^^ — )■ at T^. In 
this way 8K^'^ is not analytic at w = 0. 

Additional properties of the BCS to BEC crossover can be analyzed similarly. For example, in Appendix |B] we 
analytically evaluate the T = density structure factor in the BCS and BEC limits. Here one sees that at T = 
the density structure factor for low frequency and momentum is dominated by the gapless collective mode. As one 
crosses from BCS to BEC, this mode appears as the usual sound mode of BCS theory and evolves continuously into 
the Bogoliubov mode and eventually to the free bosonic dispersion in the BEC regime. 

One can similarly address the spin response functions following, for example the derivation in Ref. [2!]. Importantly, 
the collective modes appear only in the density response and do not couple to the spin response functions. This 
supports the discussion given in the introduction that any algebra involving both the density and spin response 
functions which decomposes these functions into separate f and J, contributions (see Eq. ^) is generally problematic, 
except in the absence of interactions. Such algebraic manipulations are not possible when the diagram sets in different 
channels are not the same. 



V. CONSISTENT LINEAR RESPONSE THEORY ABOVE T^: EXAMPLE OF PAIR CORRELATED 

STATE 

We now turn to the compressibility in a theory (of the normal phase) which includes pair correlations. It is notable 
that here too, one finds consistency with the usual Ward (and Q-limit Ward) identities, providing one restricts 
consideration to the theory originally introduced by Nozieres and Schmitt-Rink for the normal phase only. The 
NSR paper was among the first to emphasize the importance of treating pair correlations in the normal phase. Indeed 
these were discussed along with an analysis of the ground state considered here and introduced by Leggett (23 | and 
Eagles [2^. Interestingly, theories which incorporate correlated pairs which are based on this NSR scheme do not 
appear to relate to the BCS-Leggett ground state [1]. This is, in part a reflection of the rather ubiquitous first order 
transition associated with extending NSR theory below T^. Our group [131 has extensively discussed one approach 
which appears (rather uniquely) to lead to a second order transition from a different (as compared to NSR) normal 
phase into this well known ground state. However, it is more complicated than the NSR theory, and the related 
compressibility will be presented elsewhere. 

Here, in order to illustrate a fully consistent approach to linear response in the normal phase (beyond that of a 
noninteracting Fermi gas) we use the simpler NSR scheme. Even though it has been improved and reviewed many 
times (see Refs. for reviews), a full discussion on the linear response theory within NSR theory is still lacking. In 

a previous publication [sl] we have shown that by carefully choosing a set of diagrams for the vertex function, the NSR 
theory respects the usual Ward Identity associated with gauge invariance. Here we will show that the compressibility 
derived from this vertex function also satisfies the compressibility sum rule. 

We begin with a brief review of the NSR theory and its linear response theory. The self energy is S(if) — 
J2q to{Q)Go{Q - K) with to(Q) = + Xo(<3)] being the t-matrix in which xo(Q) = Ep Gq{P)Go{Q - P) is the 

pair susceptibility. The number equation is given hy n — 'i'YliK G{K) while an approximate number equation was 



> 




Figure 2: The diagrams for the vertex function of NSR theory. The first one on the right hand side is the bare vertex, the 
second one is the "MT" diagram, and the last one is the "AL" diagram. Hollow and solid dots denote full and bare vertices. 
Solid lines and wavy lines correspond to propagator of non-interacting fermions and t-matrix, respectively. Due to the two 
ways of connecting the fermion propagator inside a t-matrix, there are two AL diagrams. 

implemented in the original NSR paper Q. The response function may also be written formally as 

/^^''(^^]^, q) = 2 y r^(p + g, p)g{p + q)Y{p, p + Q)g{p) + (41) 



where the full EM vertex function must obey the Ward identity ^ so that qfj_K^'^{Q) — 0. The correction to the 
full vertex function should be consistent with that of the self energy, hence it is associated with the set of diagrams 
shown in Figl^l We have 

r^(F + Q,P)^ y\P + Q,P) + MT^(P + Q,P) + ALf'iP + Q, P), (42) 

where we identify the Maki-Thompson (MT) and Aslamazov-Larkin (AL) diagrams with 

MT'^(P + Q,P) = J2to{K)Go{K - P)-f'^{K - P,K - P - Q)Go{K - P - Q), 

K 

kU{P + Q,P) = -2 y toiK)toiK + Q)GoiK - P)GoiK - L)Go{L + Qh'^iL + Q, L)Go(L). 

L,K 

The factor 2 in the AL diagram comes from the fact that the vertex can be inserted in one of the two fermion 
propagators in the t-matrix and the minus sign comes from the fermion loop. 

To prove that the full vertex in Eq. (|42[) satisfies the Ward Identity we contract the MT and AL terms with (7^ to 
yield 

g^MT'^(P + Q,P) = Y.to(K)Go{K - P)[G^^(K - P) - G^HK - P - Q)]Go{K - P - Q), 

K 

= -[S(P)-S(P + Q)]. (43) 
(7^AL'^(P + Q, P) = -2 y to{K)to{K + Q)Go{K - P)Go{K - L)G^{L + Q)[G-\L + Q) - G-\L)]G^{L), 

L,K 

= 2[S(P)-S](P + Q)]. (44) 

In deriving the second relation, the identity Xo{K) — Xa{K + Q) = f^^^iK) — tQ^{K + Q) has been applied. One can 
then show that the Ward identity 

qt.T^{P + Q,P)^ G^\P + Q)- G^\P) + S(P) - S(P + 0) = G-i(P + 0) - G-\P) (45) 

is satisfied. Thus the linear response theory based on NSR theory with the vertex function shown in Figl2] is gauge 
invariant. 

Importantly, this same vertex also satisfies the Q-limit Ward Identity 

limrO(P + Q,P)Uo = l-^^. (46) 
By explicitly calculating the vertex function, we have 

hm rO(P + Q, P)|„^o - 1 + E *o WG^(i^ - P) - 2 V tl{K)G^{K - P)G^{K - L)Gl{L). (47) 

K L.K 



Now we evaluate the right hand side of Eq.(j46|) for NSR theory. 

, as(p) 



djjL 
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= 1 - 2 ^ tl{K)G^{K - P)Go{K ~ L)Gl{L) + ^ ta{K)Gl{K - P) 

L,K K 

= limr"(P + Q,P)U=o, (48) 

Thus the Q-Hmit WI is satisfied by the hnear response theory of the NSR theory. As a consequence, the compressibility 
sum rule is satisfied by this linear response theory. We emphasize that these two constraints (the Q-limit Ward identity 
(|46l) and the Ward identity (|45|) ') are independent constraints 0. These observations should be contrasted with the 
Hartree-Fock as well as the RPA approximations reviewed in Ref. [11 which cannot reach this level of consistency. 

Despite these findings, the compressibility is, nevertheless, ill behaved, as it diverges when the system approaches 
Tf. from above. This problem was discussed in Ref. [l^, where it was shown that the MT and AL diagrams diverge 
as T — > Tc- We note that in NSR theory, Tc is determined by the temperature where tQ^{Q = 0) = [3], as in 
the Thouless criterion in conventional superconductors (29| . In the number equation, the i-matrix resides in the 
denominator of G{K) so its divergence at Tc does not introduce difficulties. However, if one calculates dn/dfi or 
equivalently calculates the MT and AL diagrams, those expressions explicitly contain ^ in the numerator. In NSR 

theory, we have ^ = — J2q ij+^^^oil — k) + ■ ■ ■ . Due to the Thouless criterion, 1 + gx ^ o.'f for small q when 
T ^ Tc- Then the q integral behaves like J which is divergent at q — 0. 

There have been attempts to remove this divergence by including effective boson-boson interactions :30:] , but there 
has been no demonstration that gauge invariance and conservation laws are respected 31] when these repairs are 
made. One can identify more generally the problematic aspect of NSR theory. In this simplest pairing fluctuation 
scheme, where the i-matrix contains only bare Green's functions, the fermionic chemical potential is effectively the 
only parameter in the theory. Once dressed Green's functions are introduced, it becomes possible to avoid this type 
of divergence. At the same time, of course, it becomes more difficult to establish consistency with gauge invariance 
and the compressibility sum rule. Given that the compressibility sum rule is rarely satisfied (the exceptions being the 
two examples discussed in this paper) , one has a choice for how to approach a calculation of the compressibility. We 
emphasize that from the literature it appears that the more credible results for the compressibility arise via the 
thermodynamic rather than the two body correlation functions. 

VI. CONCLUSION 

Linear response theories have been an important tool for studying transport and dynamic properties of superfluid 
and related many-particle systems. The current focus in the literature on ultracold Fermi superfluids, particularly at 
unitarity, provided a primary motivation for our work which aimed to organize this subject matter and clarify the 
constraints that calibrate a linear response theory. We have seen that the challenge is to construct such theories as to 
be fully compatible with /-sum and compressibility-sum rules which reflect conservation principles. Although there 
have been some successes there are nevertheless important failures. 

In this paper we presented two nearly unique examples of fermionic superfluids which are demonstrably consistent 
with /-sum and compressibility-sum rules. We addressed these theories via the compressibility k. Important here 
was that both are compatible with the compressibility sum rule. It is useful to compare these two observations in the 
normal phase. In effect, the BCS-BEC mean field approach treated the normal phase as a normal Fermi liquid and the 
resulting compressibility is plotted in Figure [TJ Above Tc one finds very little temperature dependence. This should 
be contrasted with the behavior found in the Nozieres Schmitt-Rink approach to the normal phase, where there is a 
dramatic upturn in the compressibility with decreasing temperature. Precisely at Tc the NSR theory predicts that k 
diverges [28]. 

One can view the first of these two systems as indicating the behavior of k associated with a purely fermionic system. 
By contrast the dramatic upturn in n with decreasing T is expected for a bosonic system en route to condensation. 
Experimentally [32] the situation for unitary gases is somewhat in between these two limits. This will be an important 
topic for future research. 

Hao Guo thanks the support by National Natural Science Foundation of China (Grants No. 11204032) and 
Natural Science Foundation of Jiangsu Province, China (SBK201241926). C. C. C. acknowledges the support of the 
U.S. Department of Energy through the LANL/LDRD Program. Additional support (KL) is via NSF-MRSEC Grant 
0820054. 
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Appendix A: Detailed expressions for response functions 



The following are the EM response functions of fermionic superfluids from the CFOP theory: 

\ p p n _ f(}^^ 



+ E~ 

^ 5 [l-f{E+)-f{E-)] 



p 

A2 



iE+ + E^r^ 



E+Ep ^ co^ - {E+ - Ep)^ 



^^^^ [f{E+)-f{E-)] 



(Al) 



<3i2(w,q) = -g2i(w,q) = -iw V [(% + |^)^-^ 



V I rie- + ^P ^ 1 - /(-^P ) - /(-^P ) _ /i? _ ^P N f{E+)-f{E-) 



(A2) 



Q;3(c.,q) = = AE [(^p^ + ^p)!! - /(^p^) ' /(^p")! ^ (^^p^ - ^p-)[/(i^P^) - /(^^p")] 



p ^p 



co^-{E++Epy 



w2 - - Ep)^ 



(A3) 



Qi3(u;,q) = Qii(a;,q)=^^-^ 



- i^-)[l - f{E+) - f{E^ ^ {E++E- )[f{E+) - f{E-)] 



uj-^-{E++EpY 



w^~{E+-Ep)^ 



, (A4) 



(A5) 



<323('^>q) = -Q32i^,<l) = i 



Aw 



E^Ep 



{E++E^)[l-f{E+)-f{E-)] {E+-E^)[f{E+)-f{E^) 



LO^-{E++Ep)^ 



+ 



LO^ - {E+ - Epf 



, (A6) 



Q|3(-,q) = -Qk(-,q) = ,A V [ (^p" + ^p ' " /(^p-)i + ^^p" ' ^p )[^f p") ' /(^p-)^ 

W23^ ^32^ ^ m I w2_(^++j^-)2 ^ w^ -(£;+- i;p)2 



(A7) 



Q»(a;,q) = E [(l " 



-(1 



f+f- - A2 



+ Ep 



uj^-{E+ + Epf 



[l-f{E+)-f{E-)] 



E+Ep ^ Lo"^ - {E+ - Epf 



^^-^ [f{Et)-f{E-)\ 



(A8) 



P'P'' _ ^P + A^ 



^P^ + ^P" 



E+Ep ' uj-^ - {E+ + Epf 



[i-f{E;)-f{E-)] 



Ep Ep 



u^-{E+-Epy 



(A9) 



Q^^(-,q) = QfsM =.y^\A- ^)i_imzim^ (4 + /(^p")-/(^p) 



. (AlO) 
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Appendix B: Density Structure Factor in the BCS and BEC limits 

We first consider the BCS limit and the regime where the external frequency and momentum are small, such 
that < a; < 2 A and < (7 <C 1/^. Due to the particle-hole symmetry of strict BCS theory, Q12 = QJa = so 
K°° = + SK°°, where = Qf^ and SK°° = -Q23'932/Q22- The density structure factor is Xpp ^ XppO + &Xpp, 
where XppO = — ^Imi^Q" and 5xppO — — ^^ImtJX"'^ according to Eq. ipH]) . 

Our small frequency and small momentum limit guarantees that it is not possible to break a Cooper pair into two 
quasi-particles. Therefore XppO has no pole. Instead, Q22 determines the poles of Xpp- This leads to 

g,2(-,q) = -^(-^-c^<Z^). (Bl) 

Here -/V(O) is the density of states at the Fermi energy. The condition Q22 = yields the excitation dispersion of the 
gapless mode w = CgQ, where Cg = -7^^- Similarly, we have 

Q23(^, q) ^ ^ / d^P^ = (B2) 



p 



One then finds for the density structure factor in the BCS limit 

Xpp{uJ, q) ^Im — - — — = -^5{u: - c,9 , B3 

Cs^TT \uj — Csq + io u! + Csq + to / ZmCs 

which also satisfies the /-sum rule 

f°° 9^ 

/ rfwcJXpp(w,q) = . (B4) 
Jo ■^''^ 

Next we evaluate the density-density correlation functions in the BEC limit for different uj,q regimes. We consider 
three situations associated with (A) -j^ ^ and low frequency and momentum, (B) < 1 and low frequency and 

momentum, and (C) /2m + A respectively. Case (A) describes the deep BEC limit where a ^ and 

A/|/i| — 0. Here low momentum implies q <ti 1/C as before while low frequency means uj <C As = 2-\/|^p -t- A^, 
where As is the threshold for fermionic excitations in the BEC regime. Case (B) corresponds to a relatively shallow 
BEC regime as compared to Case (A). In Case (C), when q is sufficiently large, the system is in the very shallow 
BEC regime where /i at l/kpa — 0.553. This situation was discussed in Ref . [231 ■ The evaluation of the structure 
factor is lengthy, but straightforward. 

(A). In this case, the system is in the deep BEC limit and can be thought as a dilute gas of tightly bound molecules 
with mass ms ~ 2m. Hence, the gap A is negligible and we may approximate ~ and Qii — Q22- Here one 
finds that 



Xpp{^, q) = -5-Arcsin . " Siu - —). (B5) 

^ 1 ./-2ia-1Rl„l 



2A2(2m)l^2™^___ 2 ... . 

Note that uj = ^ = appears as an argument in the delta function, corresponding to the energy dispersion of free 
bosons with mass ms — 2rn. We note that the fermionic continuum (associated with broken pairs) does not appear 
at these low uj. 

(B). In this case, the system behaves as a weakly interacting Bose gas where the internal structure of the fermion 
pairs can not be ignored. We expand all response functions to leading order in A, q^ /2m and uj and assume _Ep ~ 
with Qii ^ Q22- The density structure factor is given by 



9 



Xpp{uo,ci)^2n^ ""-l^Siu-ujs), (B6) 



where 0;^ — \/ c^q^ + (j— )^ is the dispersion of the Bogoliugov mode. 
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(C). Since A ^ /2m + |^|, we expand the energy dispersion relation to leading order in A as = + 
A^/2^p [i^]- If, in addition, the system is in the regime < 1, we may expand the density structure factor to first 
order in A^/|/^p to obtain 



X,,(c., q) = - - - ^^)^ + --]5iu- (B7) 
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